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Abstract. We calculate the Supernova Relic Neutrino (SRN) background flux for the 
KamLAND and Super-Kamiokande (Super-K) detectors, motivated by the reduction in 
background at Super-K and new results for the star formation history (e.g., from the 
Sloan Digital Sky Survey (SDSS)). Our best estimate for the flux at Super-K is slightly 
below, but very close to the current Super-K upper limit. The Super-K upper limit is 
already inconsistent with a range of star formation histories allowed by the SDSS data. 
We estimate that the SRN background should be detected (at la) at Super-K with a total 
of about 9 years (including the existing 4 years) of data. While KamLAND is a much 
smaller detector compared to Super-K, it profits from being practically background- free 
and from its sensitivity to the lower energy supernova neutrinos. KamLAND could make 
a Icr detection of the SRN with a total of about 5 years of data. Given the small expected 
SRN event rate, we also consider the detection of the SRN in a modified Super-K detector 
with a lower threshold and reduced background where the time to detection can be reduced 
by a factor of 10 relative to the existing Super-K estimate. 
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1. Introduction 

In its death tliroes as a type II supernova (SNii) a massive star 8M0) ends its life 
emitting ~ 99% of its energy, ~ 10^'^ ergs, in neutrinos. This theoretical expectation was 
spectacularly confirmed by the detection of a handful of neutrinos from SN1987A in the 
nearby Large Magellanic Cloud jHEI- While the much larger, more heavily instrumented 
detectors such as Super-Kamiokande (Super-K) and KamLAND await the flood of neutrinos 
from another such nearby SNii, it is interesting to ask if either detector might observe the 
diffuse, isotropic flux of neutrinos from all SNii events that have ever occurred within the 
observable universe. 

The detection of this cosmological background of core-collapse supernova relic neutrinos 
(SRN) offers a new probe of SNii neutrino physics and of the high redshift Universe. The 
basic picture of core-collapse SN 11 could be tested, not only locally but also at high redshifts 
{z ^ 1). From the SRN we will also obtain estimates of the supernova rate (proportional 
to the star formation rate for M ^ 8M0) and of the metal enrichment rate which are 
completely independent of those from optical and UV surveys. Note that other tracers 
of the cosmological star formation history typically only lead to lower bounds on the star 
formation rate. Comparison with traditional methods would yield information about high 
redshift star formation complementary to that from future high redshift galaxy surveys like 
DEEP2 |S] (which targets 2^1). Here we consider the possibility of SRN detection at two 
existing neutrino detectors, Super-K and KamLAND. We also consider the detection of the 
SRN in a modified Super-K detector with lower threshold and reduced backgrounds (like 
the recently proposed SK-GADZOOKS As discussed later, detection of the SRN at 
both KamLAND and Super-K can be used to probe the cosmic star formation history at 
z ^ 1, of which little is known from traditional astronomical methods. 

The prediction of the SRN flux has been the subject of many previous investigations 
jSj. The earher study of Kaplinghat, Steigman, and Walker 2000 (KSW) reached the 
pessimistic conclusion that it would be unlikely for Super-K to detect the SRN background, a 
result driven by three factors. First there was the realization that for Super-K to detect these 
relic neutrinos, the actual flux must be close to the predicted upper bound and, in addition, 
the z/g must be nearly maximally mixed with or 1?^. However, as discussed below, recent 
results on the cosmic star formation history from observational data, including those from 
the Sloan Digital Sky Survey (SDSS), lead to a predicted SRN rate that is closer to the upper 
bound derived in KSW. Furthermore, neutrino experiments, such as SNO, KamLAND , and 
Super-K have clearly identified the large mixing angle case as the preferred solution to the 
solar neutrino problem [ZllHj, implying that the SNii Ue are (nearly) maximally mixed. The 
second factor concerned the backgrounds to the detection of the SRN. KSW estimated the 
primary background at Super-K, due to sub-Cherenkov muons (to be discussed later), by 
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extrapolating from older Kamiokande data. The estimated background to the detection of 
the SRN in the energy window (for positrons) from 19 to 35 MeV was 39 events per year. 
The latest Super-K analysis [H] reveals that the background in the same energy window is 
actually about 19 events per year. The third factor leading to the negative conclusion of 
KSW concerned the spectrum of the SRN flux. The difference between the spectra of the 
SRN and the sub-Cherenkov muon neutrino fluxes could help to separate the SRN signal 
from this background but, given the large expected background and the uncertainties in 
the supernova rate, KSW did not attempt to account for this in their analysis. However, 
considering the current status of the first two factors, it is clear that a more detailed analysis 
is now appropriate. Although with 4 years of data Super-K did not see any supernova relic 
neutrinos, they did reduce the upper limit on the SRN flux at 90% CL. to 1.2 cm~^s~^, for 
positron energies > 18 MeV jH]. This limit is more than two orders of magnitude smaller than 
the previous Kamiokande-II upper limit ^Oj and is approaching several of the theoretical 
estimates of the SRN flux. Given this very encouraging state of affairs, along with the 
advances in recent years in our understanding of the cosmic star formation history, we 
believe it is timely to attempt a realistic, quantitative estimate of the likely range for the 
SRN flux. 

For Super-K our analysis is restricted to positron energies > 18 MeV. In contrast the 
KamLAND detector, as a result of their detection technique, have much smaller backgrounds 
at energies suitable for detecting SRN and are therefore sensitive to much lower energy 
neutrinos (positrons with E > 6 MeV). Thus, even though KamLAND is a much smaller 
detector than Super-K , it profits with respect to Super-K by being sensitive to more of the 
SRN spectrum and also by its sensitivity to neutrinos from supernovae at higher redshifts 
compared to Super-K. We also present results for the SRN flux expected at KamLAND 
in the idealized limit in which KamLAND is able to utilizes its entire 1 kton volume for 
detecting SRN with 100% efficiency. We also discuss detection in the hypothetical SK- 
GADZOOKS detector, the propeties of which were described in Beacom & Vagins 
Estimates for a hypothetical detector "HyperKamiokande" fl], with a fiducial volume of 
V ~ 890 kton(«i 40Vsk), can be obtained by appropriately scaling the flux at Super-K and 
noting that the errors on the flux scale as V~^^'^. 

2. Detectors 

Super-K is a 22.5 kton fiducial mass water Cerenkov detector which can detect SRNs via 
inverse beta decay z/g + p — > e"*" + n. Below 18 MeV spallation events resulting from the 
interaction of cosmic ray muons with oxygen nuclei form the primary background at Super- 
K ^2j. Above 18 MeV there are two primary backgrounds. The first comes from low 
energy atmospheric z/g produced from the decay chain of vr^ and //^ [IBj. Below 25 MeV the 
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SRN flux from our median estimate (see below), which peaks at roughly 5MeV, exceeds 
the atmospheric Ue flux. The second background comes from atmospheric i/^ which interact 
with a nucleus to form a /x^ with kinetic energy < 53 MeV. Such fis are not detected as they 
are below the threshold for Cerenkov radiation, but their decay products are detected, 
with energy distribution given by the Michel spectrum. By performing a multicomponent 
fit to the observed Michel spectrum, as well as that for the atmospheric z/, Super-K obtained 
their impressive limit on the SRN background flux [H]. It would be ideal if there were a way 
to detect the recoil neutrons from the — p interaction (analgous to the neutral current 
detection at SNO) as was suggested in a recent proposal This neutron tagging would 
remove the and spallation backgrounds, allowing Super-K to move to a lower energy 
threshold of 10 MeV. Such a reduction in threshold will tell us more about high redshift 
supernovae and the corresponding star formation history than any other existing method. 
Regardless, for the estimates of the SRN flux at Super-K presented here we limit ourselves to 
neutrinos with energies > 19.3 MeV (positrons with Ee+ > 18 MeV). The current predictions 
from Ando, Sato, & Totani jS] show that for SK the expected flux is roughly 0.3cm~^s~^, 
which is close to the lower limit derived from Fukugita & Kawasaki Below we show how 
median estimates of the SN rate as a function of redshift can give a flux very close to the 
current Super-K upper limit. The corresponding event rate of the SRN at Super-K in this 
range is expected to be 1 to 2 per year, and a few more years of data are likely to yield a 
positive detection. 

KamLAND is a 1 kton liquid scintillation detector designed to search for evidence of 
z/e oscillations utilizing v^. from nuclear power reactors ^l^HZj. KamLAND detects v^, via 
inverse beta decay by a prompt signal from positron annihilation, followed by a ~ 200 /xs 
time-delayed neutron capture 7-ray of 2.2 MeV. The spectrum of the reactor z/g peaks 
near 3 MeV, and with neutrino oscillations is negligible above 6 MeV. The KamLAND 
collaboration analysis imposes a lower cutoff of 2.6 MeV to account for the background 
from terrestrial sources. Unlike Super-K, KamLAND can readily exclude the invisible 
muon decay background since the prompt signal from positron annihilation is not followed 
by neutron capture. The KamLAND background at energies greater than 6 MeV is from 
atmospheric z/g, which competes with our estimated SRN signal at energies ^ 25 MeV. 
KamLAND 's smaller background at low energies, in the region where the SRN spectrum 
peaks, may make it possible for KamLAND to detect the SRN background. Preliminary 
results at KamLAND for a 0.28 kton-year exposure show no z/e signal above the small 
expected backgrounds [16j. Current theoretical estimates assuming that KamLAND can 
use the entire 1 kton volume for SRN detection (Ando, Sato and Totani 2003 [S^) suggest 
that in the energy window 10 - 25 MeV, the SRN event rate is ~ 0.1 yr"^- In our 
analysis here the SRN flux for (positron) energies E > 6 MeV is considered. This lower 
energy threshold is chosen since, so far, there are no events seen above this energy at 
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KamLAND [7^. In addition, in our estimates for KamLAND we adopt an idealization in 
which KamLAND utihzes the entire 1 kton fiducial volume with 100% efficiency for SRN. 
Based on our best estimate of the supernova rate which saturates the Super-K bound, the 
expected rate in an idealized KamLAND detector is 0.4 events per year. 

3. Flux of the SRN 

The differential ffux of SRN Pg? dF/dE, depends on the magnitude and evolution (as a 
function of redshift z) of the SNii supernova rate Rsn{z), and on the Ue energy spectrum 
dN/dE, 



The energy of the detected positron, Ee+, is related to that of the Uc by Ee+ = E — 1.3 MeV. 
In the above equation the average is over the stellar initial mass function (IMF = dn^/d In M, 
where n* is the number of stars of mass M). In practice, {dN{E')/dE) is replaced by 
dN{E')/dE, with the parameters that determine the SN i/f. spectrum replaced by their 
respective IMF-averaged values. This approximation is accurate to better than 10% in the 
energy window (> 6 MeV) of interest. To compute dt/dz, we use the ACDM cosmology (i.e., 
a ffat, cosmological constant dominated model with matter density Qm = 0.3 and Hubble 
constant h = 0.7). 

3.1. Estimating the Supernova Rate 

The core-collapse SN rate at a given redshift is just the star formation rate (SFR) at that 
redshift for stellar masses larger than 8Mq. To obtain the SN rate we link the SFR to 
an observable proxy such as the UV or Ha luminosity density. The interpretation of the 
measured high redshift UV luminosity density is complicated by the fact that UV light is 
strongly absorbed by dust, while measurements of Ha, though less affected by dust, are not 
as simple to relate to the star formation rate. In addition to correcting for dust extinction, 
converting the UV light to a SFR requires correcting for the incomplete sampling of the 
luminosity function, as well as cosmological surface brightness dimming, the latter being 
true at all frequencies for any extended object [TTl UHl dl 120] 

iFiom the determination of the SFR we extract the SN rate, Rsn(^) = 
/sMo p.{z)dn,/dM{M')dM', where p*(z) is the star formation rate in average mass per 
time per comoving volume. We have limited the upper bound to 30 Mq due to the 
potential uncertainties in the neutrino ffux from SNe with progenitors more massive than 
that. Here dn^,/dlnM is the global stellar IMF, assumed constant over redshift, which is 
a good approximation provided there are no signiffcant correlations between the IMF and 
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the environment in which the stars are born. Extant evidence seem to argue against such 
correlations over the redshift range of interest (z^ 2)^]. Averaging over a Salpeter IMF 
for M > 8 M0, the SN rate is Rsn = 

assuming the star formation rate is measured 
in solar masses. The conversion factor is not sensitive to the upper limit of the average as 
long as the upper limit is larger than about 25 M0. The sensitivity to changes in the IMF 
slope is more pronounced; if we use Baldry and Glazebrook [22] (BG03; see below) best-fit 
IMF {dn^:/dln M ~ M~^'^^), the conversion factor increases by 30%. 

Following |2Sl I21j, we parametrize the SFR as p*{z) oc (1 + z)^ for 2; < 1, and 
p*{z) oc (1 + z)°' for z > 1 and use various observational proxies to estimate a, /?, and 
the normalization. From our assumption above of a redshift-independent IMF, we can then 
parametrize the SN rate as 

Rsn oc (1 + z)'^ for z <1 

oc (1 + 2)" for z>l. (2) 

In Equation 121 we have assumed that the 1 < z < 2 behavior continues to higher redshift. 

Hogg [23 has compiled measurements of the UV and Ha luminosity density, as well 
as results from measurement of the near UV emission, far-infrared and radio continuum to 
obtain the 68% c.l. limits of (3 = 2.7 ± 0.7. Results from the cosmic optical spectrum 
measurements from the Sloan Digital Sky Survey (SDSS) find limits on f3 from 2 — 3 
and a from — 2 j2^ (G03). Using the SDDS data and marginalizing over the IMF 
and evolution of the SFR, BG03 determine the local density of SFR to be in the range 
p*(0) = (0.5 — 2.9) X lO~^M0yr^^Mpc^^. The lower portion of this range for the slopes and 
normalzation is only valid, however, if the slope of the global IMF is allowed to vary near 
the BG03 best fit. This lower portion is also consistent with the results of Madau et al. 
based on the observed UV luminosity density of the whole galaxy population [2E] ■ However, 
the Madau et al. results do not take into account the correction for surface brightness 
dimming and provide a smaller average correction for dust extinction [2111, and thus permit 
the possibility that the SN rate may decline for z'^ 1. In addition, studies of extra-galactic 
background light have recently shown that current surface brightness corrections suggest 
that ~ 50% of the high redshift galaxy population is yet unobserved |27'i| . 

The choice of how to handle the dust extinction thus dominates the uncertainty budget 
for the high redshift SN rate. Although the slope of the SFR for z^ 1 is poorly constrained 
from the SDSS results, examination of high redshift UV data 20j shows that the SFR 
remains constant or even increases beyond 2; ^ 1. An increase in a would not significantly 
change the observable SRN flux at Super-K, as the energy window > 18 MeV is not very 
sensitive to the z > 1 SN rate. For example, given the parameterization in Eq. |21 we find 
that the 90 % C.L. Super-K upper bound on the flux of the relic neutrinos of 1.2 cm~^s~^ 
implies that, for P > 0.5, Rsn(O) < 2 x 10"'^ yr~^Mpc~^, independent of a. However, an 
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Figure 1. Three representative supernova rates; see eq. The sohd curve is based on 
the median results from UV luminosity density measurements as well as the SDSS cosmic 
optical spectrum. The dotted curve is an extension of this median model in which the slope 
of the z > 1 SN rate is increased from a = to a = 2; see eq. The area above the 
shaded region is ruled out by the SK flux limit. The dashed curve is a lower bound to the 
SN rate motivated by the SDSS optical spectrum. In each case a Salpeter IMF is used to 
convert from the SFR to the SN rate. 

increase in a would increase the observable flux at KamLAND (or any other relatively low 
threshold detector such as the proposed SK-GADZOOKS), due to its sensitivity to lower 
energy neutrinos. Below we discuss the extent to which a comparison of signals between 
KamLAND and Super-K can constrain the high redshift SFR. 

Motivated by the UV density studies as well as the SDSS cosmic optical spectrum, 
we choose for our "median/best estimate" model /3 = 2.5, a = 0, p^{z = 0) = 
1.6 X 10^^ M0yr~^Mpc~^. This normalization is perfectly consistent with the observations 
and is specifically chosen to saturate the Super-K result, for the adopted slopes P and a. 
This "median" SN rate (and its associated SRN background) is shown by the solid curves 
in Figures 1 and 2. Note that for this value of Rsn(2 = 0), and with an average density 
of galaxies of 0.01 Mpc~^, the present average supernova rate is expected to be ~ 1 per 65 
years per galaxy. Though the UV and SDSS analyses do permit models with higher star 
formation rates, such models are ruled out by the Super-K upper limits to the SRN flux. 
In particular, any models with values of (3 and/or Rsn(^ = 0) which exceed our median 
estimates will violate the Super-K upper bound on the SRN flux. As a specific example. 
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if the slopes were fixed at /? = 3, a = 2, and a Salpeter IMF adopted, this would require 
a SFR p^{z = 0) < 1.3 x 10"^ M0yr"^Mpc"^ in order to satisfy the Super-K flux limit. 
Thus, it is clear that the Super-K upper bound on the SRN flux already provides valuable 
constraints on the cosmic star formation history. 

It is interesting, but less certain, to estimate a reasonable lower limit to the SRN flux 
based on the above observational proxies. We choose a "lower limit" model, shown by the 
dotted curves in Figures 1 and 2, with /5 = 2, a = 0, p*(z = 0) = 0.5 x 10"^ Mgyr^-'^Mpc"^. 
Finally, in order to determine the impact on the SRN flux of increasing the high redshift 
SFR, we consider a variation to our "median/best estimate" model with (3 = 2.5 and a = 0, 
by setting a = 2 ( dashed curves in Figures 1 and 2). We note that the predictions of Ando, 
Sato, & Totani, and of Kawasaki & Fukugita [5] use SN rates which are closer to this SDSS 
inspired lower bound. 

As already mentioned, with its sensitivity to lower energy neutrinos KamLAND has 
the potential to probe the high redshift star formation history. About 40% of the flux 
in the KamLAND energy window comes from z > 1; about 10% of the flux is from SN at 
z > 2. Thus it is clear that a comparison of the fluxes at Super-K and KamLAND has the 
potential to provide valuable information about the star formation history for 1 < ^ < 2. If 
the SFR at z ^ 2 is much larger than the estimates adopted here (there are some hints in 
this direction [HIITH]), then it might be possible to probe the star formation rate at z > 2. 
Below, we quantify the information that can be gleaned about high redshift star formation 
using Super-K and KamLAND. 



3.2. SRN Spectrum 

The other key ingredient in calculating the background flux of SRN is the spectrum of 
supernova neutrinos. The supernova neutrino spectrum is typically modelled by a Fermi- 
Dirac spectrum with an effective chemical potential, rj = fiu/Ti,, 

dN_ ^ El 

dE, °^ exp{E,/T,-r]) + 1 ^ ' 

Equation ini may be used to relate the average z/g energy (-EpJ to the neutrinosphere 
temperature. Our current understanding of the SNii explosion mechanism (and the 
observations of neutrinos from SN 1987A ^1121) restricts the average energy to the range 
14 < (-EpJ ^ 17 MeV (see, e.g., |2H])- Relative to a pure Fermi-Dirac spectrum (?] = 0), for 
a distribution normalized to the total neutrino energy, positive values for rj correspond to a 
flattening of the spectrum peak, and an increase in the number of neutrinos in the tail of the 
distribution. To compromise between the extremes seen in computer simulations, we adopt 
a Fermi-Dirac distribution with zero chemical potential, characterized by (Tp^) = 5 MeV 
and (Tp^) = (Tp^) = 8 MeV. We assume that the total energy carried by each flavor of 
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neutrino is 0.5 x 10^'^ erg. The average is taken over the IMF of the SNii progenitors. This 
does not introduce additional significant uncertainties since neither the temperature of the 
neutrinosphere nor the binding energy of the neutron star depend very sensitively on the 
the mass of the progenitor and, hence, on the averaging over the IMF. The most recent 
simulations of Raffelt et al. jSHI including muon and tau transport models suggest that the 
IMF- averaged energies may not be hierarchical, but that (-Ep^) ~ (Ep^^). For a spectrum 
normalized to the same total energy output in all neutrino species, such differences do not 
significantly change the results for the fluxes. 

3.3. Neutrino Mixing 

Neutrino oscillations are of direct relevance to the spectrum of background SRN. If z/^ and/or 
Uj- mix with the spectrum of the resulting z/g will be harder. Such higher energy z/e are 
easier to detect. From the analyses of solar neutrino and KamLAND data, mixing between 
the t'e(i^e) and I'^ip^) flavor eigenstates is favored, with 0.27 < tan^6'0 < 0.94 at 3cr jS]- 
In terms of mass eigenstates, the MSW solution implies AttIq = m\ — m\ > The 
atmospheric neutrino data are consistent with near-maximal vacuum mixing between 
the v^iv^) and frii^r) flavor eigenstates: sin^(26') > 0.82 at 90% C.L. However, unlike the 
solar neutrino and KamLAND data, the atmospheric neutrino data do not constrain the 
sign of the corresponding mass-splitting, Am^^^ = |m| — m^. Here, two mass hierarchies 
are possible; a normal hierarchy, with > m2 > rui, or an inverted hierarchy, with 
m2 > mi > 777,3. As will be noted next, the choice of hierarchy is important for any neutrino 
mixing inside the supernovae. 

Due to the high density inside supernovae, there is the possiblity that before reaching 
the surface of the supernovae, the neutrinos may have been mixed by the MSW effect |33j . 
In supernovae, there are two MSW resonant density layers, corresponding to the solar and 
atmospheric mass splittings. Whether the resonant density occurs in the v 01 v channels 
depends on the mass hierarchy j31|. For a normal hierarchy, both resonant layers are in the 
V channel. However, the flux that arrives at the surface of the supernova is an incoherent 
mixture of mass eigenstates. These mass eigenstates travel separately to the surface of the 
earth and, in the absence of earth matter effects, the final z/g flux will have a fraction sin^ 6q of 
the original z/^ flux. In contrast, for an inverted hierarchy, the higher density resonant layer 
is in the z/ channel. Deep inside the supernova the z/g flavor corresponds to the lowest mass 
eigenstate, whereas in vacuum Vr is the lowest mass eigenstate. For adiabatic propagation 
inside the supernova, and for |Ue3| ~ (but not exactly zero) as imphed by the reactor 
neutrino data jlH], Ve remains the lowest mass eigenstate as it exits the supernova, while 
Vj. remains the heavy mass eigenstate. The result is that for an inverted hierarchy the z/g 
observed on Earth were all "born" as Vr- 
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4. Results 

Having assembled the necessary ingredients, we are now in a position to calculate the SRN 
flux at Earth and to estimate the event rates at KamLAND and Super-K. For our standard 
cases we assume a normal mass hierarchy (no mixing inside the supernovae) and choose 
tan^^0 = 0.46 (sin^(^0) = 0.31), in the middle of the solar neutrino and KamLAND range 
(see ^3.3|) . For Super-K we find 

0.3 < F{E > 18 MeV) < 1.2 cm'^ s'^ (4) 

while for KamLAND, 

1.9 < F{E > 6 MeV) < 7.8 cm'^ g-^. (5) 

The upper and lower limits correspond to the solid and dashed curves in Figure |21 
respectively. For an inverted mass hierarchy, the detectable flux for E > 18 MeV increases 
by ~ 50%, as this interval samples the high-energy tail of the SN neutrino spectrum 
(locally and, especially, at higher redshifts). If this inverted hierarchy were realized, a lower 
supernova rate would be required in order to remain consistent with the Super-K upper limit. 
An inverted hierarchy is less important for the interval E > 6 MeV, which is sensitive to 
more of the (zero-redshift) spectrum, and for our SN neutrino spectrum and SFR parameters 
there is a negligible change in the detectable rate. Note that the median/best estimate SFR 
(solid curve in Figure saturates the current Super-K upper bound of 1.2 Ue cm~^s~^ [Oj, 
suggesting that with a further reduction in the background, Super-K may detect the SRN 
background (see below for further discussion). 

In order to calculate event rates corresponding to these neutrino fluxes, the relevant 
cross sections and energy thresholds are required. Neglecting the smearing of events due 
to finite energy resolution, the event rate for positrons produced by the SRN with neutrino 
energies between Ei and E2 (assuming 100% efficiency) is: 

R(Ei,i?2) = Np/ a{E)—dE, (6) 

JEi+l.SMeV dE 

where E is the Ue energy, Np is the number of protons in the fiducial mass of the detector and 
cr(i?) is the Ue cross section on protons Here we take E2 = 82 MeV. In the left panel 
of Figure ^ we show the results for the integral event rates for KamLAND, Super-K, and 
SK-GADZOOKS. The results for Super-K cover the range of SN rates allowed by combining 
the SDSS and SK limits, while for KamLAND we show the integral event rate for our best 
estimate. We estimate a Super-K event rate 

1.1 < Rate {El > 18 MeV) < 3.6 events/yr, (7) 

and a KamLAND event rate 

0.1 < Rate [El > 6 MeV) < 0.4 events/yr. (8) 
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Although KamLAND is a much smaller detector, it does have the advantage of 
no background and sensitivity to lower energies, permitting it to sample higher redshift 
supernovae and more of the energy spectrum of lower redshift supernovae. This effect 
is illustrated in Figure © where the fractions of the SRN fluxes contributed by SNii at 
different redshifts are compared for KamLAND {Ei = 6 MeV) and Super-K {Ei = 18 MeV). 
However, if the SFR has been underestimated at higher redshifts, KamLAND will detect 
even more relic neutrinos. For example, increasing the slope of the high-z supernova rate 
from a = to 2 results in a ~ 36% increase in flux for Ei > 6 MeV, but only a ~ 4% change 
in the event rate for the Super-K window Ei > 18 MeV. 




Figure 2. Left: Estimates of tiie supernova relic integral flux for positron energies 
E > £^„iin, as a function of i?min- As in Figure 1, the solid curve is our "median SDSS" 
model with a normalization chosen to saturate the Super-K upper bound. The dotted 
curve, also saturating the Super-K upper bound, shows the effect of increasing the slope 
of the z > 1 SN rate. The Super-K upper bound is indicated by the cross. The dashed 
curve is from an estimated lower bound to the SN rate based on the SDSS results. The 
Kamland and Super-K threshholds are shown by the vertical lines. Right: Estimates of the 
SRN Dc integral event rates for positron energies E > ii^min, as a function of -Emin. The 
dashed, dotted, and upper solid curves correspond to those for Super-K in Figure ^ The 
lower solid curve corresponds to the solid curve in Figure ^ for the integral event rate at 
KamLAND. Crosses mark the levels of the predicted event rates for our best estimate for 
each detector. 
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Figure 3. The fraction of the rehc neutrmo flux contributed by SNii in each redshift 
interval. The dotted histogram is for Super-K, while the solid histogram is for KamLAND. 



What if both KamLAND and Super-K should detect the SRN? This would open up 
the possibility of learning about the star formation history for z > 1, which is currently 
not well understood. A differential comparison of the KamLAND and Super-K signals, the 
latter of which comes primarily from z < 1 (Fig. (jH))), would contain information about the 
z > 1 SRN rate and, therefore, about the z > 1 supernova rate. It is informative to define 
the quantity 

r = - 1, (9) 

where Rk,sk and Vk^sk are the event rates and fiducial volumes at KamLAND (K) and 
Super-K. The parameter r faithfully captures the effect of changes in Rsn(2;) on the ratio of 
the flux in the 6-18 MeV window compared to that in the > 18 MeV window. Note that 
r does not depend on the normalization Rsn(O); it is only sensitive to the shape of the Rsn 
- redshift relation Rsn (2;). Herein lies the sensitivity to supernovae at 2; > 1. If the slope 
of the z > 1 SFR is increased from a = to 2, r changes by roughly 25%. For example, 
for the ranges of event rates discussed above, 1.8 < r < 2.2. However, changes in a are, 
to some extent, degenerate with changes in /?. The advantages for studying the evolution 
of the cosmic SFR provided by the detection of the SRN at KamLAND will come from 
combining such information with other, more traditional, astronomical observations of the 
cosmic SFR which can further restrict the range of f3. Such analyses will provide unique 
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information about the SFR at high redshift. 



5. Prospects For Detection of the SRN 

If the SRN background is indeed at the level of our median estimate (middle curve of Figure 
12} , how long will it take for Super-K and KamLAND to detect the SRN? To answer this 
question we must estimate the size of the error bars on the flux of SRN obtainable at Super- 
K and KamLAND, as a function of the number of years of data accumulated, T. This 
expected error can then be compared to the likely ranges for the SRN-expected flux. Our 
model for the error only includes the statistical errors, which in fact is a good approximation 
jH]. The data, given the parameters which determine the SRN and the atmospheric neutrino 
fluxes, is expected to be Poisson distributed. We neglect the correlations between bins due 
to energy resolution effects. If the bin sizes are larger than about twice the energy resolution, 
the correlations will change our results at the ~ 10% level. 

The recent Super-K analysis jH] lists the detection efficiency for the SRN events as 
47% for El < 34 MeV and 79% for Ei > 34 MeV after full data reduction. We use an 
efficiency of unity for Kamland. We modify the event rate for the SRN and the backgrounds 
in accordance with this prescription for efficiency. The event rate in each energy bin is 
denoted by Sai where i refers to the energy bin. a identifies whether the event rate is due to 
the SRN, invisible muon decay, or the atmospheric neutrinos; /„ is the corresponding flux 
normalization. Assume that each bin accumulates Ni events (after full data reduction) in 
T years of runtime. The Poisson measurement error is a"^ = Ni. We define 

F(a,/?) = rS,^^^. (10) 

The error on /„ is then given by [F~^{a,a)Y^'^. Assuming the model described above is a 
good fit to data (as demonstrated by the Super-K results IH]), we forecast the errors on fa 
by replacing N^ by its expected mean value, TTi^Sai- 

With the above efficiency corrections, for the Super-K window Ei > 18 MeV the 
expected 1-cr error on the SRN event rate with T years of data is ~ 11/ VT events/yr. We 
note that with a total of 8 — 9 years of data, Super-K will be able to make a 1-a detection 
if the SRN fiux is close to our estimate from median SFR parameters (solid line in Figure 
12)). For KamLAND , the primary background for the SRN is from atmospheric i^e(^'e), 
which has negligible expected rate in the range ^ 25 MeV. To determine this background 
at KamLAND, we correct the results from Super-K for the number of target protons in 
the 1 kton KamLAND volume. We determine that KamLAND can make a 1-a detection 
of the SRN fiux with ~ 5 kton-years of data, if the fiux is close to the solid curve in Figure 

m 
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These estimated times to detection are long because the expected signal is close to 
the background when the detector is big and high threshold (Super-K) or small and low 
background (KamLAND). We note that if Super-K can sufficiently reduce the background 
from spallation and sub-Cerenkov /i^ decays, it can become an optimal SRN detector, 
opening up the E > 10 MeV window for detection of the SRN j4ij. Assuming 80% efficiency 
for detecting the SRN for E > 10 MeV (M. Vagins, private communication), and a factor 
of 5 reduction in the backgrounds 4J, we determine that Super-K will make an immediate 
detection of the SRN background, at minimum a 1-a detection with less than one year 
of data. This opens up the possibility that we could constrain the space of fundamental 
SN II parameters which affect the spectrum of neutrino flux from a SN ii significantly. The 
general strategy would be to relax some assumption (motivated by simulations) about the 
neutrino spectrum and test them with SRN data. 

6. Conclusions 

Stimulated by the Super-K collaboration's P] remarkable reduction in their backgrounds 
for detection of the SN relic neutrinos [E > 18 MeV), by the absence of significant 
backgrounds to the SRN detection at KamLAND |2|(for E > 6 MeV), and by recent 
progress in pinning down the cosmic star formation history from new observations including 
the Sloan Digital Sky Survey, we have calculated the expected fluxes and event rates at 
Super-K and KamLAND. The SRN flux presented here, based on the SN rate derived from 
the star formation data, saturates the current Super-K upper bound, suggesting that Super- 
K may be close to observing the SRN. For our median model (see Fig. |21), we conclude 
that the predicted flux at Super-K saturates the current experimental upper limit, with 
a corresponding event rate of 3.6 yr^^. In this case Super-K will be able to make a 1-a 
detection by doubling their data set, i.e. with a total of about 9 years (including the existing 
4 years of data). If Super-K does see the SRN background flux, assuming that KamLAND is 
100% efficient and can use the entire 1 kton volume for detection of the SRN, the associated 
event rate expected at KamLAND is small, ~ 0.4 yr^^. However, in this ideal scenario, the 
lack of any significant competing backgrounds will permit a 1-a detection of the SRN with 
only 5 years of data. 
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